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We study the production of sterile neutrinos in the region T ~ Mw in an extension beyond the 
standard model with the see-saw mass matrix originating in Yukawa couplings to Higgs-like scalars 
with masses and vev's of the order of the electroweak scale. Sterile neutrinos are produced by the 
decay of scalars and standard model vector bosons. We obtain the index of refraction, dispersion 
relations, mixing angles in the medium and production rates including those for right-handed sterile 
neutrinos, from the standard model and beyond the standard model self-energies. For 1 < Mw /T < 
' 3 we find narrow MSW resonances with k < T for both left and right handed neutrinos even in 

I absence of a lepton asymmetry in the (active) neutrino sector, as well as very low energy {k/T <C 

^ . 1^1) narrow MSW resonances in the presence of a lepton asymmetry consistent with the bounds 

O ' from WMAP and BBN. For small vacuum mixing angle, consistent with observational bounds, the 

j /' j ' absorptive part of the self-energies lead to a strong damping regime very near the resonances resulting 

, in the exact degeneracy of the propagating modes with a concomitant breakdown of adiabaticity. 

We argue that cosmological expansion sweeps through the resonances, resonant and non-resonant 
sterile neutrino production results in a highly non-thermal distribution function enhanced at small 
momentum k < T, with potentially important consequences for their free streaming length and 
f~ l ■ transfer function at small scales. 
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^ . I. INTRODUCTION 

> , 

OO . In the concordance ACDM standard cosmological model, the Universe today is composed approximately by 70% 
of a dark energy component responsible for the acceleration, about 25% of dark matter (DM) and about 5% of 
ordinary matter (baryons). In this scenario the (DM) component is cold and coUisionless, and structure formation 
proceeds in a hierarchical "bottom-up" manner: small scales become non-linear and collapse first and their merger 
and accretion lead to structures on larger scales This is a consequence of the fact that cold dark matter (CDM) 
particles feature negligible velocity dispersion leading to a power spectrum that favors small scales. In this hierarchical 
scenario dense clumps that survive the merger process form satellite galaxies. Numerical simulations of structure 
formation with (CDM) predict many orders of magnitude more (DM) sub-haloes than observed low luminosity dwarf 
. ^ galaxiesB [1, 3, [p , H- These simulations also yield a density profile that increases monotonically towards the 
K> , center [sTItI. la. lol. [ToT flit p{r) ^ r"'' , 7 = 1 corresponds to the Navarro-Frenk- White (NFW) profile, but steeper 
5^ ' profiles with 7 ~ 1.2 have been found recently in numerical simulations [3| ■ These density profiles accurately describe 
] clusters of galaxies but there has been recent observational evidence that seems to indicate a shallow cored profile 
instead of cusps in dwarf spheroidal galaxies (dShps) which are deemed to be (DM) dominated [T^ ITsl IT3. ITsl ITa. Iit\ . 
This core vs. cusp controversy is still being debated, and recent arguments suggest that the interpretation of the data 
is subject to assumptions and modelling J_S] . Recently yet another discrepancy between the predictions of ACDM and 
observations has been revealed, the "emptiness of voids" , possibl y re lated to the overabundance problem p^. 

Warm dark matter particles (WDM) were invoked [20I [21I [2^. |23|| as possible solutions to the core vs. cusps and 
the overabundance problems in satellite galaxies. (WDM) particles feature a non-vanishing velocity dispersion with 
a range in between (CDM) and hot dark matter (HDM) leading to a free streaming scale that cutoffs power at small 
scales thereby smoothing out small scale structure. If the free streaming scale of the (WDM) particles is smaller than 
the scale of galaxy clusters, the large scale structure properties are indistinguishable from those of (CDM), but may 
affect structure at small scales [l^l, thereby providing an explanation of the smoother inner profiles and the fewer 
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satellites. A small scale cutoff in the (DM) power spectrum may also explain the apparent smallness of galaxies at 
z - 3 found in ref.[25l. 

Although the interpretation of cores in (dSphs) may be challenged by alternative explanations, and the missing 
satellite problem could be resolved by astrophysical mechanisms such as complex "gastrophysics" , and recent simula- 
tions suggest that the dynamics of subhalos is not too different in (WDM) and (CDM) models^], there is an intrinsic 
interest in studying alternatives to the standard (CDM) paradigm. 

Any particle physics explanation of (DM) involves extensions beyond the Standard Model (SM), allowing quite 
generally, both (CDM) and (WDM) candidates. 

Sterile neutrinos, namely SU{2) singlets, with masses in the ~ keV range may be suitable (WDM) candidates [13, [H, 
[Hj [13, HH [Hj [Hj Si mi and may provide possible solutions to other astrophysical problems [3l|. The main property 
that is relevant for structure formation of any (DM) candidate is its distribution function after decoupling [3J, |36[ 
which depends on the production mechanism and the quantum kinetic evolution from production to decoupling. 
Sterile neutrinos may be produced by various different mechanisms [13, [H, [IS [13, [U, [H, [H|, among them non- 
resonant mixing, or Dodelson-Widrow (DW)[23, [H, [l^l has been invoked often. However, there seems to be some 
tension between the X-rayfsS] and the Lyman-a forest datajs^, [i^l leading to the suggestion [4ll| that (DW)-produced 
sterile neutrinos cannot be the dominant (WDM) component. 

A phenomenologically appealing extension of the (SM) with only one scale has been recently proposedfsol [ssl [svj . 
In this model sterile neutrinos may be produced by the decay of a gauge-singlet scalar with a mass of the order of the 
electroweak scalej31^.32„^, 42l_In this scenario sterile neutrinos are produced and decouple at a temperature of the 
order of the mass or the scalar [321 [ssl. l43j . 

Recentlyji^] the quantum kinetics of production and decoupling of ^ keV sterile neutrinos in these models was 
studied with the result that production via the decay of the gauge-singlet scalar leads to a non-thermal distribution 
function that favors small momentum. This result was combined with an analytic method to obtain the transfer 
function during matter domination recently introduced in ref. [4^ . This method reveals the influence of the distribution 
function of the decoupled particle upon the power spectrum and free streaming lengthji^. The results of re f.[43| p oint 
out that sterile neutrinos produced via the decay of gauge-singlet scalars in the model advocated in refs. jSOl. l3ll. [3^. 
[ssl . [4^ . [4^ 1 yield smaller free streaming lengths and an enhancement of power at small scales as compared to those 
produced by the (DW) mechanism. Combining the results for the distribution function of sterile neutrinos produced 
via scalar decay with abundance and phase space constraints from dwarf spheroidal galaxies (dShps)[33] yields a 
narrow window for the mass of sterile neutrinos [43|: 0.56 keV < Ms < 1.33 keV. The robustness of this bound has 
been confirmed in ref. [i^l , but there may be some tension with recent analysis of the Lyman-a forest with non-thermal 
populations [4^, although the results in this reference relate mainly to resonant production. 

Recent observations of the X-ray spectra from the Ursa Minor dwarf spheroidal galaxy with the Suzako satellite [isj 
suggest that sterile neutrinos with masses in the keV range with mixing angles 9 ~ 10~^ remain viable candidates as 
main dark matter constituents, a result that seems to be confirmed by those of ref. [13]. 

Motivation and objectives: The clustering properties of dark matter candidates depend on the free streaming 
length which determines the scale below which power is suppressed. When the (DM) particle of mass Ms has become 
non-relativistic, the free streaming length is approximately given by 



A 



fs 



where p is the (DM) density and the average is with the distribution function of the decoupled (DM) particle. 
Distribution functions that favor small momenta lead to smaller free streaming lengths and more power a small 
scales m [3 [H. 

The study in ref. 14311 revealed that the non-resonant mechanism of sterile neutrino production by scalar decay 
advocated in ref. 13, 12, [H, [13 leads to a non-thermal distribution function that favors small momenta with impor- 



tant consequences for structure formation and remarkable differences with sterile neutrinos produced by the (DW) 
mechanism [23|, whose distribution function is that of a thermal relic that decoupled while relativistic, but multiplied 
by an overall factor [23|. This overall factor in the (DW) distribution function only affects the abundance, but for a 
fixed (DM) density p the resulting free streaming length is that of a neutrino of mass Ms decoupled at the sterile 
neutrino decoupling temperature. For a fixed mass and relic density the non-thermal distribution function from the 
production mechanism studied in ref. (3^ . [33l . [43^ yields a smaller X fs and more power at small scales than in the 
(DW) mechanism without modifying the large scale power spectrum. 

In the extension beyond the standard model (bsm) advocated in ref .lsol . [H, [H, \^ \^ sterile neutrinos mix with 
active neutrinos via a Yukawa coupling to the standard model Higgs[30j whose expectation value yields a see-saw 
mass matrix. The diagonalization of this see-saw mass matrix yields interaction vertices between vector bosons and 
the sterile-like neutrino. This is important because the distribution functions being a function of the energy, are 
necessarily associated with mass or energy eigenstates, not flavor eigenstates. 
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The study in ref . [H, reveals that sterile neutrinos are produced and decouple at a temperature of the order 

of the mass of the scalar, which in the model of refs.fsO, [H, [ssl. l42j| is of the order of the Higgs mass. 

At this temperature the charged and neutral vector bosons are present in the medium with large abundance, 
comparable to that of the scalar. Their decay into the sterile-like neutrinos will therefore contribute to their total 
abundance and distribution function. This is one of the main observations in this article. The coupling of the charged 
and neutral vector bosons to the sterilc-likc neutrino is suppressed by the (small) mixing angle, but since the standard 
model couplings are much larger than the Yukawa couplings of the scalar to the sterile neutrino, the question is 
whether the decay of vector bosons may lead to a substantial contribution to the production rate of the sterile-like 
neutrinos. For Ms ^ keV and the expectation value of the Higgs-like scalar in the range ~ 100 GeV the Yukawa 
coupling Y ~ 10~^, the production rate via this process oc Y^, whereas the contribution from Z,W decay would 
be expected to be c>c atuSin^(6'), with 9 the mixing angle. For 9 ^ 10~^^^ the production rate of sterile-like 
neutrinos via vector boson decay can be of the same order of or larger than that from scalar decay. This observation 
suggests that sterile neutrino production via the decay of vector bosons in the medium may be competitive with the 
production via scalar decay. 

At high temperature and or density the mixing angle is modified by medium effects [i^, HO, Ell, [13, [13, Issl ]. 
therefore the first step towards understanding whether vector boson decay contributes to the production of sterile-like 
neutrinos is to obtain the in-medium correction to the mixing angles. 

A more general aspect of sterile neutrino production via vector boson decay at T ~ A1\y is that both the index of 
refraction (real part of the self-energy) and the production rate determined by the absorptive part (imaginary part 
of the self-energy) are of 0{Gf)- This is in contrast to the usual situation at temperatures much smaller than the 
electroweak scale when the index of refraction is of 0{Gf), but the absorptive part is of 0(0"^). 

Althou gh t he finite temperature and density corrections to the index of refraction have been obtained for T <C 
Mm/ z (53I l54i [55} . to the best of our knowledge the study of the self-energy, the index of refraction (real part) and 
absorptive part (width) at T ~ Miy,z has not been carried out. 

To be sure, upon the diagonalization of the mass matrix, standard model interaction vertices with the sterile-like 
neutrino lead to production processes via both charged and neutral current interactions such as II —^T'iV2 ; // — > 'viV2 
with charged leptons (1) or quarks (f) and vi ^ Va',V2 ^ i^s (active and sterile respectively) for small mixing angle. 
These processes are of 0{a^ sm^ 9) and while they will eventually become important for T <C Mw,z when the 
population of vector bosons in the medium becomes <^ a^, these are formally subleading in the weak coupling at 
T ~ Mw,z- 

Therefore at T Mw^z vector boson decay is the leading production mechanism from weak interactions. 
Our objective is to provide a comprehensive assessment of sterile neutrinos as potential (DM) candidates imple- 
menting the following program: 

• Obtain the production rates and mixing angle in the medium from the quantum field theory model at T ^ Mw^ 
studying the possibility of MSW resonances to determine whether sterile neutrinos are produced resonantly or 
non-resonantly. 

• Obtain and solve the kinetic equations describing the production and decoupling of sterile neutrinos using the 
rates and mixing angles obtained from previous step. 

• The asymptotic solution of the kinetic equation yields the distribution function after freeze out, which determines 
the abundance and the free streaming length. This distribution function is input in the program described in 
refs.ji^ll^] to obtain the transfer function and power spectrum. 

In this article we carry out the first step of this program. We implement methods of field theory at finite temperature 
and density developed in refs.[5^. [sgI . IstI [ssI . [53 | to obtain the mixing angles in the medium and production rates 
both from scalar and vector boson decay. 

Results: We study a simple extension of the standard model with one active and one sterile neutrino to extract 
the robust features in a simpler setting. Both active and sterile neutrinos are considered to be Dirac, this is to include 
the possibility of a lepton asymmetry hidden in the (active) neutrino sector (Majorana neutrinos cannot be assigned 
a chemical potential), and to allow us to study the production of left and right-handed neutrinos. 

We obtain the dispersion relations, index of refraction, mixing angles and production rates in the medium from the 
self-energy contributions from standard model (sm) and beyond the standard model (bsm) interactions. The see-saw 
mass matrix that mixes them emerges from the Yukawa couplings to Higgs-like scalars with masses of the order of 
Mw,z that acquire expectation values also of this order. We focus on the temperature region T ~ Mw where vector 
and scalar bosons are present in the medium with large thermal populations. The decay of both the scalar and vector 
bosons contribute to the production of sterile neutrinos. Our main results are: 
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• We find one MSW resonance even in absence of a lepton asymmetry. For 1 < Mw/T < 3 this resonance 
is in the low momentum region 0.2 < k/T < 1 and well within the regime of validity of the perturbative 
expansion. Including a lepton asymmetry in the active neutrino sector consistent with the data from Wilkinson 
Microwave Anisotropy Probe (WMAP) and Big Bang Nucleosynthesis (BBN)[62'|, we find two low energy MSW 
resonances, the lowest one is a consequence of the lepton asymmetry that occurs at k/T <C with |^| being the 
lepton asymmetry parameter. In the region of interest for this study for small vacuum mixing angle consistent 
with the observational bounds from X-ray datafi^ these resonances are very narrow. We find resonances also 
for positive energy, positive helicity, namely nearly right-handed neutrinos. 

• At the resonances the propagating frequencies become exactly degenerate in striking contrast with the quantum 
mechanics of neutrino mixing wherein there is level repulsion at the resonance. This exact degeneracy at the 
resonance entails the breakdown of adiabaticity. It is a distinct consequence of the absorptive part of the 
self-energy and leads to a strong damping regime. 

• The form of the standard model contribution to the production rate is similar to that from scalar decay found in 
ref. [i^ . We argue that cosmological expansion will lead to a rapid crossing of the narrow resonances resulting in 
both resonant and non-resonant sterile neutrino production. In particular nearly right handed sterile neutrinos 
are produced by the decay of Z°, vector bosons. Their distribution functions after freeze-out will be highly 
non-thermal with a distinct enhancement at small momentum k < T and perturbatively small population. This 
low momentum enhancement of the non-thermal distribution function is expected [43l. [43 | to have important 
consequences: a shortening of the free-streaming length (smaller velocity dispersion) and an increase of the 
transfer function and power spectrum at small scales. 

• We find a consistent range of parameters for which there is a resonance for positive helicity, positive energy 
neutrinos, namely nearly right-handed at T ~ Mw- The general field theory framework allows a systematic 
study of the properties for both helicity states, including the helicity dependence of mixing angles and production 
rates. 



II. THE MODEL 



The extension of the standard model presented in ref. P, '33', "425 generalizes the proposal of the z^-MSM of ref. [30|, |37[ 
and is also a generalization of the model presented in ref. [60]. These models include three SU(2) singlet (sterile) 
neutrinos which couple to the active neutrinos via a see-saw mass matrix. The generalization of ref. [H, [H, H^] gives 
a mass to the sterile neutrino via a Yukawa coupling to a Higgs-like scalar field which could be the neutral Higgs 
component, or another scalar whose expectation value is of the same order as that of the (sm) Higgs boson, therefore 
this type of extension features only one scale. 

We study a simplified version of these models by considering only one sterile and one active neutrino. In the 
usual see-saw mechanism an off-diagonal Dirac mass matrix for the active species is considered along with a diagonal 
Majorana mass for the sterile neutrino jsd IstI fsol. fs lll52l| . However, instead of considering a Majorana sterile neutrino, 
we allow for Dirac mass terms for all species. This generalization allows to study simultaneously the possibility of a 
lepton asymmetry in the (active) neutrino sector for which a Dirac field is required, along with the possibility of a 
right-handed component leading to potentially relevant degrees of freedom within the same simple model. 

Our goal is to extract generic and robust features of the production rates and mixing angles in the medium 
along with a reliable estimate of sterile production rates. The generalization to three species can be done relatively 
straightforwardly (but for the complications associated with dealing with larger mixing matrices), and the case of a 
Majorana neutrino is regained straightforwardly by projection. 

We consider a model with one active (va) and one sterile (vs) (an SU{2) singlet) Dirac neutrinos, described by the 
Lagrangian density 



£ = CsM ^v,i^Vs-Yx VsHU - Y2 Vs<^iys + C[<i>] + h.c , (H.l) 

where 

^ • (n.2) 

/ is the charged lepton associated with and H'^, are the components of the standard model Higgs doublet, and 
is a real scalar singlet field whose expectation value gives a Dirac mass to the sterile neutrino. The Lagrangian 
density describes the kinetic and potential terms of <I>. 
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In unitary gauge we write 



(11.3) 



and consistently with the single scale assumption of the i^-MSM: ~ ($) are of the same order of magnitude (the 
weak scale ) and that their masses are also of the same scale. In fact our analysis is quite general, and this assumption 
will only be invoked for a quantitative assessment. The Lagrangian density pi.l[l becomes 



-SM 



Ma/3 l^fJ - Yi VsCrVa - Y2 Vs^Vg + £[($) + ip] + Yi.c \ a, (5 ■ 



a, s 



where 



(11.4) 



M = 



m 
m AL 



(11.5) 



Introducing the "flavor" doublet {va, Vs) the diagonalization of the mass term M is achieved by a unitary transforma- 
tion to the mass basis (z^i, 1^2), namely 



where 



In the mass basis 



Mm = U-^{0) ; MU{e) 



{''j)=u{e)[ll) ; u{e)=( '^"'(^^ 



cos(26i) 



Ml 
M, 



M, 



Ml 



sm{9) cos{e) ) ' 



2m 



Ms - {Ml + 4m2] ■ 



(II.6) 



(II.7) 



(II., 



We focus on a see-saw with Mg 

Ml - 



keV ^ TO therefore 
,2 



TO 



— ; M2~Ms 



M2 



< 1 . 



(II.9) 



Taking {H'^) ~ ($) the small mixing angle entails that Yi <C 1^2 which results in self-energy corrections from the a 
exchange are subleading as compared to those from the ip exchange. For example taking ($) ^ and for a ^ keV 

sterile neutrino it follows that 



Y2^10-^:^Yi ; sin(20) - Fi/Fa . 



(11.10) 



However, we can alternatively consider a pre-determined see-saw mass matrix and set Yi = Y2 = which corre- 
sponding to a simpler extension of the standard model that posits a mass matrix that originates beyond the standard 
model. 

Our goal is to obtain the dynamical aspects of sterile neutrinos in the medium, mixing angles, dispersion relations 
and damping rates, which determine the production rates. These are obtained directly from the solution of the 
equations of motion including the self-energy corrections in the medium. The one-loop self-energies require the 
neutrino propagators in the medium in the mass basis, since the mass eigenstates are the true propagating states. For 
6 <^ 1 the mass eigenstates J^i ~ I'a ; 2^2 ~ , and the active neutrino reaches equilibrium at T > 1 MeV via the weak 
interactions, whereas the sterile neutrinos are not expected to equilibrate. 

This argument, however, hinges on the smallness of the vacuum mixing angle, but in a medium the mixing angle can 
become very large, and if there are MSW resonances the roles of the medium eigenstates may be reversed. Whether 
there are MSW resonances and the medium mixing angle becomes large can only be answered a posteriori. 

Therefore we assume that the mass eigenstate vi is active-like, and features a Fermi-Dirac distribution function, 
whereas for z/2 the propagators are the vacuum ones. Furthermore, it is possible that if there is a large lepton 
asymmetry it may be stored in the neutrino sector, whereas the asymmetry in the charged leptons equals the baryon 
asymmetry and can be neglected. Hence the Fermi-Dirac distribution functions in the vi propagator includes a 
chemical potential. 

In our study we explicitly separate the fermionic and bosonic contributions to the self-energies to assess the consis- 
tency of the assumption that the eigenstate "1" is active-like. 
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III. EQUATIONS OF MOTION 

The effective Dirac equation in the medium is derived with the methods of non-equihbrium quantum field theory 
described in [s^ls^lsoj. We follow the approach presented in refs. [sgI. [STj and introduce an external Grassmann- valued 
source that couples linearly to the neutrino field via the Lagrangian density 

whence the total lagrangian density is given by £ + £5. The external source induces an expectation value for the 
neutrino field which obeys the eff'ective equation of motion with self-energy corrections from the medium [59|. 

The equation of motion is derived by shifting the field i^^ — 4'a + ^ V'q = Wa) imposing (^^) = order by 
order in the perturbation theory [sl, [sj, [s^ . Since the self-energy corrections to the equations of motion require the 
neutrino propagators, we obtain the equation of motion in the mass basis. 

Implementing this program up to one loop order, we find the following equation of motion for the doublet in the 
mass basis -0"^ = (^1 1 "^2), it is given by 



where I is the identity matrix, = diag{Mi, M2) is the mass matrix in the mass basis, L = (1 — 7^)/2 is the left- 
handed chiral projection operator, S*^ is the (local) tadpole contribution from the (sm) neutral current interaction, 
(see Fig. ([T])) .YZ!^{x — x',t — t') and SJ|^(af — x',t — t') are respectively the real-time retarded self-energies from 
(sm) and (bsm) (scalar) interactions. Introducing the space-time Fourier transform in a spatial volume V 

^{x, t)^^Y.j ^^0^*"^ e-^'°*^(fco, k) (ni.3) 

A: 

and similarly for the self-energy kernels and the source term, the equation of motion in the mass basis becomes 



(70^0 - 7 ■ - + Y.'Z L + S,,„(fco, k)L + S6.™(fco, fc ) 



Vi(fco,fc) --rXfco,fc) . (III.4) 



The space-time Fourier transform of the retarded self-energies (not the tadpole) feature a dispersive representation 

nk..k)^'-r duJ^^^th%. (III.5) 

A. One-Loop Self-Energy 

We focus on the temperature region Mz,w,iy,tp ^ in which using the unperturbed thermal propagators for the 
scalar and vector bosons is valid|611- In section (|V D[) wc show that perturbation theory is valid for k ^.otw T for 
Mw,a,Lp ^ furthermore for k <^ Mw our results reproduce those found in the literature for T <C [ssl. [s^ and 
the perturbative expansion is reliable for Mw > 2T. 

The (sm) charged and neutral current contributions to the self-energy in the mass basis are depicted in Fig.((T]). The 
latin indices i, /c = 1, 2 refer to the mass basis fields and the label / in the intermediate fermion propagator in the 
charged current diagram in Fig.([T|) refers to the charged lepton associated with the active neutrino. The contributions 
from scalar exchange (bsm) in the mass basis are depicted in Fig. ([2]). 

(SM) neutral currents: The tadpole contribution in the mass basis is given by 

= u-\e){l l]u{e) , (iii.e) 



where 



5^ f d^l r - \ _ 9^ 



^1 

'7t2 



e = f. (111.7) 



^ This expression corrects a typographic error in ref. |56|). 
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i j i k j 




i f j 

FIG. 1: Standard model contributions to the self-energy "Esm- The indices i, k,j = 1,2 corresponding to mass eigenstates, the 
index / for the intermediate fermion line in the charged-current self-energy refers to the charged lepton associated with the 
active neutrino. 



, I , 1 ^ 

i k j 



FIG. 2: Beyond the standard model contributions to the self-energy T.bsm- The indices i,k,j = 1,2 corresponding to mass 
eigenstates. The dashed line is a scalar propagator either for a ox ip 



In this expression n^,ni, are the Fernii-Dirac distribution functions for neutrinos and antineutrinos respectively, and 
we have neglected the contribution from the asymmetry of the charged lepton and quark sectors since these are 
proportional to the (negligible) baryon asymmetry. We allow for a lepton asymmetry stored in the neutrino sector. 
A recent analysis [62i| from the latest WMAP and BBN data suggests that |^| < 10"^. 

The neutral current diagrams that contribute to the one-loop self energy feature two different terms corresponding 
to the intermediate neutrino line being either vi or i'2- As argued above, for small mixing angles i>i ~ i>a and weak 
interactions equilibrate these mass eigenstates with the medium, therefore their finite temperature propagator features 
the Fermi-Dirac distribution function (with a chemical potential allowing for a lepton asymmetry). However, 1/2 ^ Vs 
will not equilibrate with the medium since their coupling to the environmental degrees of freedom is suppressed by 
at least two powers of the (small) mixing angle, therefore 1/2 features a vacuum propagator. The one loop diagrams 
are shown in Fig. ([3]) where the superscripts (1) and (2) are used to specify the intermediate neutrino propagator vi 
and V2 respectively. 



z^ 




i 1 j i 2 j 



y.(2) 

FIG. 3: Neutral currents contribution to the one-loop retarded self-energy Esm. The indices i,j = 1,2 and the indices 1,2 
denote the corresponding mass eigenstate in the intermediate state. 
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g... (7.. 




^(1) ' ' y(2) 

FIG. 4: Scalar exchange contributions to the one-loop self-energy Tibsm- The indices i,j = 1,2 and the indices 1, 2 denote the 
corresponding mass eigenstate in the intermediate state. 



In the mass basis we find for the neutral current contributions shown in Fig. ([3]) 



S„c(fco, k) 



'(0)I]«(fco,A;) +sin2(0) T.^}{ko,k)\ U-\9) 



1 




uio) 



(IIU 



(sm) charged currents: the charged current one-loop self energy is shown in Fig. ([T]), since the intermediate 
state is a charged lepton we find in the mass basis 



'Ecc{ko,k) ^ 'Ecc,sm{ko,k)U ^{9) 



1 




uio) 



(III.9) 



where Scc,sm(fcoi k) is the usual standard model one- loop self-energy in thermal equilibrium. 

(bsm) scalar exchange: The scalar exchange contributions to the self-energy are shown in Fig. ([3]). For sin^(6') <C 1 
we find 



-'bsm 



(fco, k) 



(0) ^i'Hko,k) + sin2(^?)sW(fco, k) + cos2(0)E(2) (fc^^ fc) 



\0)^<^^Hko,k) u-\e) 



1 




u-\e) 




1 



u{e) 



U{0) 



(III.IO) 



Summary of self-energies in the flavor basis: The structure of the self-energies (to leading order in sin^ {6 j) 
(|III.6miI.10[) indicates that they are diagonal in the flavor basis. In this basis the total self-energy is given by 



Saa(fco,fc) 

T,ss{ko,k) 



(III.ll) 



where 



Saa(fco,fc) = [s*-Kcos2(0)I]W(fcQ^fc)+sin2((?)E(f2(fco,fc) + S™(fco,fc)]i + cos2(^^ , (in.12) 

^ss{ko,k) = cos2(0)S^i)(fco,fc)+sin2(0)EW(fco,fc)+cos2(^?)s(,2)(fco,fc) . (11113) 

Since in the (sm) contributions we have explicitly factored out the left-handed projector L, the remainder contribu- 
tions to the (sm) self-energies 'Snccc are those of a vector-like theory. The (bsm) contributions feature both chiralities 
since we have considered a Dirac mass term for the sterile neutrino, a left-handed Majorana mass term can be obtained 
by neglecting the right-handed contribution. We consider the regime T 3> Mi 2, ruf and kQ,k ^ Mi_2, ™/, where to/ 
stand for the charged lepton masses, therefore we can safely neglect the mass terms and consider the propagators of 
massless fermionic fields. 
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In this regime the general form of the (sm) self-energies with vector boson exchange, either charged or neutral 
currents is written in dispersive form as in eqn. piI.Sp withjSl, Ist} 



/mSsm(a;, k) = ng 



2 

sm 



i2n) 



(III. 14) 

where F stands for the fermionic species in the intermediate state. For vi and charged lepton n^? is the Fermi-Dirac 
distribution function, whereas for 1/2 it is np ~ since the "sterile" neutrino does not thcrmalizc with the medium. 
For the (bsm) contributions, the general form for scalar exchange is 



/ml]6sm(tj, fc) = ttF^ 



d^ 



where 



dpodqo5{uj - po - qo) I - npipo) + Nsiqa) PF{po,P)PB{qa,q) 



CC 

^ NC 



(III.15) 



(III. 16) 



^ 2cos(0^) 

and Y — Yi,Y2 for a and exchange respectively. The spectral densities are respectively (for massless fermions) 

Pf{P(),P) = \{^°^^'~) (^(Po -p) + ^(7° + 7 • p S{Po+p), 
PB{qo,q) 

The projection operator 



2W, 



\S{qo-Wq)-6{qo + Wg)] ; Wg = Vq^ + AP . 



(III.17) 
(III.18) 



and 



npiPo) 



Pi^y{qo,q) 



1 

5(Po-A')/T -I- 1 



; q"^ 



{q\q) 



; "-f(po) = 1 - ny{-po) ■,NBiqo) 



eQo/T _ I 



(III. 19) 



(III.20) 



We have allowed a chemical potential for the neutrinos (only for vi ^ Va) to include the possibility of a lepton 
asymmetry in the (active) neutrino sector. 

In the expressions above, the masses for the scalars or vector bosons are M^^ip, Mz,w as appropriate for each 
contribution. All the self-energies share the general form 



E(fco, k) = 7° A{kQ,k) - 7 • k B{kQ,k) , 



(III.21) 



the detailed expressions for the imaginary parts of the (sm) and (bsm) contributions are given in the appendices. 

In particular, for the neutral current tadpole B{kQ,k) — and A(kQ,k) can be recognized from eqn. (jlll.Tp . 
Combining pil.ll[) with this form we write the self-energy matrix in the flavor basis as 



{ko,k) L + 'Sbsmikojk) 



7° AL(/co,A:)-7-kBL(fco,A:) 



L- 



7° Ai?(fco,A:)-7-kBfl(fco,fc) 



In the flavor basis these matrices are of the form 



A{k„,k) = 



AaaikQ,k) 







Ass{ko,k) 



; B(fco,A:) = 



Baaiko,k) 







B,,{ko,k) 



R. (III.22) 



(III.23) 



where the matrix elements are obtained from the expressions (jIII.12IIII.13p . 

The equations of motion for the left (L) and right (R) handed components are obtained by multiplying the equation 
of motion pil.4p on the left by the projectors R and L respectively. 
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It proves convenient at this stage to separate the Dirac spinors into the left ■>pL and right ipfi, handed components 
and to expand them into hehcity eigenstates[5^, namely 



h=±l 



and 



h _ 



(III.24) 



(III.25) 



where the left Lp and right C, handed doublets are written in the flavor basis, and are eigenstates of the helicity 
operator 



1 

1 



(III.26) 



namely, 



(III.27) 



To leading order in weak and Yukawa couplings, and neglecting a commutator [M, S] because it is higher order in 
these couplings, we find in the flavor basis for both the left and right-handed component doublets 



[kl ~ k^)! + (fco - hk) (Al + KEl) + (fco + hk) {Ar - KEr) 



It 



(III.28) 



where M is the mass matrix in the flavor basis and the inhomogeneities in these equations are obtained by projection 
and using the corresponding equations, we need not specify them as they are no longer used in our study. 

In absence of interactions, for the left-handed component a positive energy solution corresponds to /i = — 1 and a 
negative energy solution to h = +1 with the opposite assignment for the right-handed component. 

In the flavor basis 



2 \ sm(2t') cos(20) 



(III.29) 



where 



M 



■ {Ml + M|) ; 5M^ = M| - 



and Mi,2 are given by eqn. (jll.Sp . 

It proves convenient to define the combinations 



Sh{ko, k) 



(fco + hk) \{Ar - /iBfi)^^ + {Ar - KUr)^ 

(ko - hk) \{Al + hBL)^^ + {Al + HBl)^ 



and 



Ahiko,k) 



(fco + hk) 



(fco — hk) 



5A'P 



{AR-hBR)^^-{AR-hBR) 

{Al + h^L) {Al + hBL) 



where we have suppressed the arguments. The equation of motion pil.28p can now be written as 



(III.30) 



(III.31) 



(III.32) 



(III.33) 
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where the inverse propagator is given by 



where 



and 



-Chiko,k) Dhiko,k) 
Dh{ko,k) Ch{ko,k) 



Ph{ko, k) 



(cos(26l) + Ah{kQ, k)^ + sin2(20) 



Chiko^k) 



Dh{kQ,k) = 



(cos(20) + A^(fco,A:)) 



Ph{ko, k) 



sin(26') 

Ph{ka,k) 



(III.34) 



(III. 35) 



(III.36) 
(III.37) 



We note that if Afi{kQ,k) were real, then Ch{kQ, k) = cos,{26'^-^{kQ, k)) and Dhikg, k) = s\ii{29^{kQ,k)) with 9^{kt^,k) 
the mixing angle in the medium for the different helicity projections and as a function of frequency and momentum. 



B. Propagator: complex poles and propagating modes in the medium 
From (|III.34P we read off the propagator projected onto helicity eigenstates 

(E,i(fco,fc) = 

where 



2(a;j(A:o, k) - /3/i(fco, k) ) 2[ah{kQ, k) + (3h{ko,k) 



(III.38) 



Th{ka,k) 



-Ch{ka,k) Dh{ko,k) 
Dh{kQ,k) Ch{ko,k) 



(III.39) 



I 2 

ah{ko,k) = kl- k^ + -Shiko,k)- M , 
Ph{ko,k) - ^5M^ph{ko,k) . 



(III.40) 
(III.41) 



If Ah{ko, k) given by eqn. (jIII.32|) were real, the propagator pil.38|l would be diagonalized by the unitary transfor- 
mation 



u^{et{ko,k)) 



cosi0t{ko,k)) sm{0'^{ko,k)) 
-sin(C(fco,fc)) cos(C(fco,fc)) 



(III.42) 



leading to 



1 



u-\em)G{ko,k)u{e„ 



a{ko, k) + (3{ko,k) 




(III.43) 



i{ko,k) - /3{ko,k) / 



where we have suppressed the helicity argument for simplicity. However, because Afi{ko,k) features an imaginary 
part determined by the absorptive part of the self-energies, there is no unitary transformation that diagonalizes the 
propagator. However, since the imaginary part is perturbatively small the expression (jIII.43p clearly indicates that 
the pole for a = /? corresponds to the mass eigenstate 2, namely a sterile-like neutrino state, and the pole for a — ~(3 
corresponds to the mass eigenstate 1, namely an active-like state. 
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We note that in absence of interactions, namely Sh = ; Ah = it follows that 

a + 13 = k^-k^- , 
a- 13 = kl-k^ - M| 



(III.44) 
(III.45) 



The propagating eigenstates in the medium are determined by the (complex) poles of the propagator pil.38p , which 
again correspond to ah(kQ, k) = ±j3h(k(j, k). 

Before we analyze the complex poles, it proves convenient to separate the real and imaginary parts of a, (3. For this 
purpose and to simplify notation, we suppress the label h and the arguments /cq: k in these quantities, and we write 



S = SR, + iSi ; A = An. + iAj 



(III.46) 



where the subscripts i?, / stand for real and imaginary parts respectively. Furthermore, we define the mixing angles 
in the medium solely in terms of the real parts of the self-energy (index of refraction), namely 



, cos(26') + A7? . 

cos(26'„) = ; sin(26l„ 

Po 



sm{29) 

PO 



where 



Pa 



cos(26') 



sin2(26') 



An MSW resonance occurs whenever cos(26'm) = [43, [SO, [mI, [S^] , namely when 

Afl = - cos(26') . 



(III.47) 



(III.48) 



(III.49) 



We emphasize that the mixing angle in the medium 9rn and po depend on helicity, fco, k. In terms of these quantities 
we find 



P = ~~^Po r [ cos((/)) + i sin(0)] = (3^ + iPi , 



where 



(l-7')' + (27COs(2f?„))' 



; 7 



and 



sign(7COs(26lm))< -arctg 



27COs(26'„ 



1 — 7' 



TT 1 



0(1 - r) + 2 - 2arctg 



27COs(2&„.) 



1 — 7 



e(f -1) 



(III.50) 



(III.51) 



(III.52) 



This form is similar to that obtained in a model of oscillations and damping with mixed neutrinos studied in ref. ^63| , 
and suggests two distinct situations: a weak damping case for It"! < 1 and a strong damping case for I7I > 1. These 
cases will be analyzed below. 

Zeroes of a + (3: We are concerned with the ultrarelativistic limit k ^ M| ^ Mf. Just as in the usual 
case it is convenient to introduce the average or reference frequency 



uj{k) ^yk^ + M . 

The poles are near Lj{k), therefore write 

ko = uj{k) + {ko — cJ(fc)) 
keeping only the linear term in (/cq — w(fc)) we find 



a + /3 - 2uj{k) 



ko-ni{k) + iTi{k) 



(III.53) 
(III.54) 

(III.55) 



13 



with 



ni{k) = uJ{k) - — 

ri(fc) 



Sr + SM^por cos((?!)) 



4cZJ(fc) 

Si + SM'^porsm{(t)) 



- kQ—uj{k) 



4cZJ(fc) 

Zeroes of a — (3: Proceeding in the same manner, we find 



(III. 56) 
(III.57) 



ko—ij(k) 



with 



(3 - 2U{k) 



n2ik) = LJ{k) - — 



4u{k) 



ko - n2ik) + ir2{k) 



Sr - SM'^por cos{(t>) 



T^ik) 



Auj{k) 



Si - SM^por sin(0) 



(III. 



(III.59) 



(III.60) 



=IU(fe) 



From (|III.55IIII.58|) it is clear that the propagator in the medium features two Breit-Wigner complex poles corre- 
sponding to the two propagating modes in the medium. 

In the expressions above we have only focused on the positive energy modes. The expressions for the negative 
energy modes may be obtained from the following relations which are consequences of the imaginary parts of the 
self-energies and the dispersive representation valid both for scalar and vector boson exchange (jIII.5|) , 

lmA{-kQ,k; fi) = ImA(fco, fc; -/i) ; RcA{-kQ,k; fi) = -ReA{ko,k;~p,) , (III. 61) 

lmB{-ko,k;p.) = -lmB{ko, k; -p) ; ReB(-fco, fc; ^) = ReB(fco, A:; -^) . (III.62) 

These properties can be read-off the explicit expressions for the imaginary parts of the self-energies given in the 
appendix equations (jA.mA.3|) for the standard model contributions and equations (|B.mB.3| for the scalar exchange 
contributions. The matrices A are extracted from the coefficient of 7" and B from the coefficients of 7 • fc in the 
self-energies respectively. The relations for the real parts follow from the dispersive representation piI.Sp . 
In what follows we use the ultrarelativistic approximation 

Lij{k) ~ k 



M 

~2k ' 

In the limit of interest k/T < 1 with Mi < M2 Ms ^ 0{keV), the region k <T 
wide window in which the ultrarelativistic approximation is reliable. 

We note that the difference in the real part of the pole position in the ultrarelativistic limit becomes 

6M^ 



^2{k) - ni{k) 



2k 



-porcos{<p) . 



(III.63) 

0(100 GeF) corresponds to a 



(III.64) 



From the expression (jIII.52p for I7I > 1 it follows that when an MSW resonance occurs, namely for 9^ — 7r/4 
resulting in cos(0) = and the real part of the poles become degenerate. This is in striking contrast with the 
quantum mechanical description of mixed neutrinos where no level crossing (or complete degeneracy) can occur. 
Indeed the degeneracy is a consequence of the fact that the self-energy is complex and only occurs when damping is 
strong in the sense that I7I > 1. 

The degeneracy near an MSW resonance for strong damping will necessarily result in a breakdown of adiabaticity 
during cosmological evolution. We analyze below the conditions required for this phenomenon to occur. 

Furthermore, as discussed in refs.|32|, [4^ decoupling and freeze-out of sterile neutrinos of neutrinos produced via 
scalar decay occurs near the electroweak scale, and it will be seen consistently that vector boson decay yields a 
production rate with a similar structure as for scalar decay therefore a similar range of temperatures in which sterile 
neutrino production by this mechanism is effective. 

Perturbation theory is reliable when the change in the dispersion relations (positions of the poles in the propagators) 
is small. In the relativistic limit the (bare) poles correspond to ko = k (for positive energy particles), therefore 
perturbation theory is valid for k 3> {^1.2 — k) ; ri^2 namely k ^ 5](fc, k) where E is any of the self-energies. In 
the next section we obtain explicitly the self energies and in section (jVPp we assess the regime of validity of the 
perturbative expansion. 
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C. Helicity dependence: right-handed sterile neutrinos and standard model interactions 



We have purposely kept the general form of the self-energies and propagators in terms of the helicity projections 
h = ±1. In the non-interacting massless case, positive energy left-handed particles correspond to h = —1 and negative 
energy left-handed correspond to /i = 1, with the opposite assignment for right-handed particles. For the massive but 
ultrarelativistic case the mass term yields corrections to the handedness-helicity assignment of 0(M^/fc^). 

h = 1: Neglecting subleading terms of 0{M jk"^) that multiply (bsm) right-handed contributions in the ultra- 
relativistic limit, we obtain 



S{k) 
A{k) 



2k 



2k - 

Jap . 



(Az 



- (Al 



0. 



(III.65) 
(III.66) 



h = 1: In this case the corrections of 0{M /fc^) multiply (sm) left-handed contributions, which may be of the same 
order of the (bsm) right-handed contributions. We find, 



S{k) 

m 



= 2k 



2k r 



A, 



(Az 



^(Az^+Bz.)^^ 



(III.67) 
(III.68) 



The terms proportional to M /4fc^ only receive contribution from the standard model self-energies, whereas the 
right-handed components only originate in the contributions beyond the standard model which are suppressed by 
much smaller Yukawa couplings. However the last contribution in ()III.68P from (sm) interactions may he of the same 
order as the (bsm) contributions for a relevant range of k. To see this note that Azj,Bz;, ~ Y2 ~ 10~^^, whereas 
Ai, Bz, ~ (7^ ~ 0.4 therefore with M ^ KeV and k < 100 GeV, it is clear that both contributions (bsm) and (sm) are 
of the same order. 

The point of maintaining the helicity dependence throughout is that for the case of sterile neutrinos, namely the 
propagating modes "2" in the medium, the exchange of standard model vector bosons yields a contribution to the 
positive helicity and positive energy components, namely the right-handed component, which could be of the same 
order of the (bsm) contributions for small k which is a region of interest for sterile neutrino production. 



IV. REAL PARTS: MIXING ANGLES AND MSW RESONANCES 

The dispersion relations (real parts of the poles) and the mixing angles in the medium are determined by the real 
parts of the self-energy, namely the "index of refraction" . Whereas the neutral current tadpole contribution pil.7p is 
real, the real part of the other contributions is obtained from the dispersive form pil.sp . namely 



ReE(fco.fc) - / dujvi 



Iml](ti), k ] 
LU — ko 



(IV.l) 



In general the real part must be obtained numerically and is a function of three parameters /cq, fc,/z which makes its 
exploration a daunting numerical task. However, progress can be made by focusing on the "on-shell" contribution, 
namely setting /cq ~ fc, and neglecting the dependence on /i, which is warranted in the whole region of k, T of interest, 
but for k/T, \^J-\/T <C M/T, in which case we provide below an accurate approximate form. 

In obtaining the real parts we consider only the finite temperature contribution, because the zero temperature part 
is absorbed in the renormalization of the parameters in the Lagrangian. 

Scalars (bsm): For the real part of the scalar (bsm) self-energy we find for ko — k ; /i = 



ReY,bsm{k, k) = 



16^' 



k M\ ^/k M 



7 • k 



k M 
T' Y 



Bb 



k M 
T'Y 



(IV.2) 
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FIG. 5: The functions Af{k/T; M/T); Ab{k/T- M/T); Bf{k/T; M/T); Bh{k/T; M/T) as a function of k/T for M/T = 1, 2, 3 

for ReEi 



where Af; Bf and Ab; Bb are the fermionic and bosonic contributions respectively and Y — Y\,2 for cr, iy9 exchange. 
Figs. © show A/; 5/ and Ab; Bb for A//T = 1, 2, 3 as a function of k/T. 

For Sct,^ the intermediate fermion Hne corresponds to a sterile-Hke neutrino, therefore for these contributions we 
must set Aj = 0; Bf = 0, under the assumption that the sterile neutrino population can be neglected and the 
propagator for the internal line is the vacuum one. For the mixing angle the relevant contribution is A — B. Figs. 
display Af - Bf and Ab - Bb for M/T = 1, 2, 3 as a function of k/T. 





FIG. 6: The functions Af{k/T;M/T)-Bf{k/T;M/T);Ah[k/T;M/T)~Bh{k/T;M/T) as a function of fc/T for M/T = 1,2,3 

for the case of scalars (bsm). 



We note that the fermionic and bosonic contributions A/, Ab are qualitatively very similar and the same property 
holds for Bf, Bb. Therefore neglecting the fermionic contributions both for E^^^^ does not affect the results and the 
conclusions in a substantial manner. 
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This observation confirms that the general results presented below are robust even when the neutrinos "1" are not 
thermalized and their propagators are the vacuum ones. 

Although an analytic form for the full range of fco;fc;/i is not available, we obtain an analytic expression for the 
relevant case fco/T, k/T, fi/T ^ M/T ^ 1. We find to leading order in the smaU ratios ko/T; k/T; ^ = ^/T, 



12 



Tw^koT^ 
120 M2 



1 + F[M„/T] 



'j-k 



7n^ kT' 



■2 r 



360 M2 



1 + J[M^/T] 



(IV.3) 



Y2rp2 

ReE<,')(t„,4) = -^. 



12 I TT' 



77r2 ko T- 



2 r 



120 M2 



1+F[M^/T] 



-j-k 



360 M2 



1 + J[M^/T] 



(IV.4) 



V'2t2 



77r2 fcp 
120 M2 



F[M„/T] -i-k 



7^2 fcr2 

360 M2 



J[Ma/T] 



(IV.5) 



y27"2 

ReE(,2)(fco,A;) = -2^. 



7^2 kor^ 



F[MJT] 



120 M2 ^ 



7 • fc 



77r2fcr2 



360 M2 " ^^•^'^ 



(IV.6) 



where 



120 r°° fl2 



120 /"^ o2 



. 2 



(IV.7) 



These functions are displayed in Fig. they are 0(1) in the region of interest Ma-^ip ^ T. 
A comprehensive numerical study of Af, Ab, Bf, Bb confirms the validity of the above approximations for fco 
fc,^ = for fc/T< 1. 

Vector bosons (sm): Similarly, for the real part of the (sm) self-energy we find for fco = fc ; ^ — 

— 7 • k 



167r2 



, , k M 
Bb[ -: — 

T T 



(IV.8) 



where we use the same definition, namely Af] Bf and Ab: Bb arc the fermionic and bosonic contributions respectively. 
Figs. (0) show Af] Bf and Ab] Bb and Fig.® shows Af{k/T] M/T)-Bf{k/T] M/T)]Ab{k/T] M/T)-Bf{k/T] M/T) 
for Af/T = 1, 2, 3 as a function of k/T. 

Just as in the (bsm) case analyzed above, we note that the fermionic and bosonic contributions Af] Ab are qualita- 
tively similar and the same holds for Bf] Bb. Again this observation confirms that our results are robust, independently 
of whether any of the neutrino modes is thermalized. 

We also obtain the analytic forms for ReY,sm{ko] k) for ko/T, k/T,fi/T ^ M\y,z/T ^ 1. To leading order in these 
small ratios we find 



ReS«(fco,A:) 



4M2, 



^ 1 + ^ 



7^2 fcoT2 

60 Ml 



l + G[Mz/T] 



'7-fc 



77r2 kT 
180 M| 



2 r 



l + G[Mz/T] 



(IV.9) 



ReE(il)(fco,fc) 



9^ 

4M2, 



77r2 kpT^ 
60 Ml 



G[Mz/T] 



7 • /c 



7 2 L.rp2 

180 Ml ^ J 



(IV. 10) 



ReS^ 



liko.k) 



2M2, 



7^2 kpT^ 
60 M2, 



1 + G[Mw/T] 



"f ■ k 



7^2 fcT2 

180 M^ 



1 -1- G[M„./r] 
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FIG. 7: The functions Af(k/T; M/T); Ab{k/T; M/T); Bf{k/T; M/T); Bb(k/T; M/T) as a function of k/T for M/T = 1, 2, 3 

for the standard model contributions (sm) with /i = 0. 
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FIG. 8: The functions Af{k/T; M/T) - B f{k/T; M/T); Ab{k/T; M/T) ~ Bb{k/T; M/T) as a function of fc/T for M/T = 1,2,3 

for the standard model contributions (sm) with fi = 0. 



In the charged current contribution we have neglected the asymmetry of the charged lepton because it is of the 
order of the baryon asymmetry. In the above expressions 



G\m] 



120 

7^ Jo 



dq ^Nb{W,) 



« 4 



NsiWg) 



1 



W„ 



(IV. 12) 



This function is depicted in Fig.®, it is 0(1) in the region of interest T ~ Mz,w- 

The validity of these approximations for k{) = k, ji = Q is confirmed by the numerical analysis of Af, Ab^ Bf, Bb for 
/c/T<l. 

It is remarkable that the leading order in k^/T, k/T but for Mw,z ~ T reproduce the results of references [ssl 
which were obtained in the low energy limit T,fi <^ Mw,z- The numerical analysis carried out for ko = k; ii = 



18 




confirms that for M /T 3> 1 the range of validity of the lowest order approximation in k/T increases and merges with 
the results given above in eqns. (|IV.3IIIV.TT|) up to k/T ^ 1. 



A. Mixing angles and MSW resonances: 



As shown in detail in the previous sections, the mixing angle in the medium 9^ determined by the relations (jIII.47[) 
depends on fcoi ^ and the helicity h. On the mass shell of the propagating modes in the medium we can replace ki^ ^k 
in the expressions for the real part of the matrices ReA ; RcB for A/j, namely the real part of eqns. (jIII.66IIII.68[) 
for h — ^1 respectively. For fi — 0; ko — k and general k, M the fermionic and bosonic contributions to the real parts 
of the (bsm) self-energies are given by eqn. (IIV.2|) where the fermionic (A/, Bf) and bosonic {Ah, Bh) contributions 
are depicted in figs. The real parts of the (sm) self-energies are given by eqn. pV.8|l and the fermionic and 

bosonic contributions depicted in figs.([7]l8]). 

These figures distinctly show that the contributions Af, Ah and Bf, Bh for (bsm) and (sm) self-energies are quali- 
tatively the same, with only a quantitative difference in the amplitudes. A remarkable result is that these functions 
change sign. In particular the combinations Af — Bf, Ah — Bb which enter in change sign at a value of k/T that 
depends on the ratio M/T. For M/T 1 these differences vanish at k/T ~ 0.2. A numerical exploration reveals that 
the sign change persists until M/T ~ 3 but occurs at monotonically larger values of k/T. This behavior is shown in 
the figures above. We find that for M/T > 3 the change in sign occurs for k >> T or does not occur at all. On the 
mass shell ko ^ k and for fi = this study reveals that Aji is negative in a wide region of momentum for M/T < 1. 
This fact entails that there are MSW resonances near the momentum regions where the coefficient functions change 
sign, even in absence of a lepton asymmetry. To understand this important point more clearly let us study the case 
h = ^1 separately. 

h — — 1: In this case A is given by eqn. (jIII.66p . furthermore from eqn. (jIIL13|) it follows that (Al — ]Bl)ss is 
determined by the (bsm) contributions which are suppressed by small Yukawa couplings Y < 10^^ as compared to the 
(sm) contributions. Therefore the (bsm) contribution can be neglected and A^ is determined by the (sm) contributions 
given by eqns. pV.8p . furthermore approximating cos(6') ~ l;sin(6') ^ in eqn. (IIII.12P and (5M^ ~ M^, we find (for 




Af 



1 



2 cos{9m) 

Taking as representative T ~ 100 GeV; M. 



k Mu 



T' T 



Abi 



(k Mw\ Mw\ 



k Mw' 



k M 



w 



T' T 



Ab(K'-^)-Bf(K^)-Bb(K^ 
\T T J \T T J \T T 



keV it follows that 

2 



167r2 



T 
M, 



— -2.7x10". 



(IV. 13) 



(IV. 14) 



Figures ^ show that for M-i^^z/T < 3 there is a region in k/T in which the bracket in pV.13|) is negative and there 
is a value {k/T)c that increases with M/T at which the bracket vanishes, for example from the Fig. ^ we find 
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(fc/T)c 0.2; 0.45; 1 for M/T 1, 2, 3 respectively. For k/T < {k/T)c the bracket is positive (for fi = 0) whereas for 
k/T > {k/T)c it is negative, therefore there is a value of k/T at which the resonance condition (|III.49[) is fulfilled. 
Since the coefficient of the bracket is « 10^^ (eqn. (jfV.14ll ) and the terms inside the bracket are of 0(1) for k/T < 1, 
and cos(6') ~ 1 it follows that the MSW resonance occurs for a value of k/T such that the bracket ~ 10~^^ namely 
for k/T ~ {k/T)c- The large coefficient (|IV.14I) results in a very narrow MSW resonance as can be seen as follows, 
expanding A^j near [k/T)c as 



+ 



ft > , 



(IV. 15) 



where k > 10^^ for M-^^z/T < 3 (see Fig. (|5])) and approximating cos(20) ~ 1 we find 



sin2(26l„0 



(*)-(*).-0 



e = sin(26')/K . 



(IV. 16) 



For example taking sm{26) ~ 10^^ [48] it follows that e < 10^-'^^ which makes the resonance very narrow. During 
cosmological expansion the ratio M/T{t) increases with the scale factor, while the ratio k/T (with k the physical 
momentum) is fixed. Therefore, for a fixed value of k/T < 1 as M/T increases the resonance is crossed very sharply. 

h = 1: To assess the possibility of MSW resonances for ft, = 1 we need the real part of (|III.68|) . From (|III.12|) 
and pV.2p it follows that (A/j — JiRjaa oc ; (A/j — Bfl)^^ oc Y2 , since I2 ^ Yi we can neglect the first term 
(corresponding to cr-exchange). Similarly in the term (A^ — Bi^)^^ we neglect the contribution from (T-exchange and 
approximate cos{9) ^ 1 ; sin(0) ~ in piI.13p . hence only E^^"* contributes to E^s- Furthermore, approximating 
5M^ ~ M ^ ~ A/2 we finally find for ft = 1 ; ^ = 0, 




(IV. 17) 



where in the first line the Ab\ Bh are (bsm) displayed in figs. ([5]). 

We note that with T - lOOGeV, ~ KeV the value of the Y2 (see eqn. (|IlT01) is such that Y2T/MS - 0(1), 
therefore fig ^ (right panel) suggests that the (bsm) contribution may yield an MSW resonance in the region 
k/T ^ 0.15 ; ^ T, where the (bsm) contribution Ab — Bh is positive and large. Since g'^ /12d>-n'^ ~ 3.4 x 10"'' and 
Af + Bf;Ab + Bb ^ 0(1) for k/T < 1 it follows that the (sm) contribution to is subleading for k/T < 1 and 
the (bsm) contribution may lead to an MSW resonance in this region depending on the parameters of the extension 
(bsm). 

7^ ; k/T < M/T - 1 : 

The above results are valid for fi — 0, for fi ^ a, full numerical evaluation of the real parts of the kernel is not 
available, however, the bounds on the lepton asymmetry from ref.fi^] suggest that l/x/Tj < 0.02 ^ 1 and we can 
obtain a reliable understanding of the influence of the lepton asymmetry (in the neutrino sector) by focusing on the 
region of fc/T < 1, in which we can use the results pV.3llV.6|) for (bsm) and pil.7p along with pV.9IIIV.ll"|) for (sm) 
and approximate cos(6') l;sin(6') ~ in piI.12IIII.13|) and 5M^ ~ M^. 

For ft = — 1 again we neglect the (bsm) contributions to Aii(k) in pil.66p . and for fJ,/T; k/T ^ 1 we obtain. 
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and the resonance condition (|III.49p can be fulfilled for ^ > when the bracket in (jIV.18|) approximately vanishes, 
namely for 



(I) 



25 ^ 
56^ 



k - 0.05 fi , 



(IV.20) 



where we have used G{AIw,z /T) ^ 1 for T ^ Mw, a result that can be gleaned from Fig. For ^ > t his MSW 
resonance occurs for antineutrinos (namely fco = —k), a result that follows from the relations (|III.61|III.62p . 



Similarly, for ft, = 1 and /i/T, k/T <C 1, we obtain 
A«(^) - -(^) 
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Obviously, there is a competition between (sm) and (bsm) contributions in eqn. pV.2ip . When T 
J(1),F(1),G(1) ~ 1 and {Y^T'^)/M'^ ~ 1. Therefore, the (bsm) contribution to Afl,(fc) is 
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and the (sm) contribution to Afl(fc) reads 



^ \k) 192 VT/ 240 
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The resonance happens for Afl(A:) — 1, namely 
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Obviously, one is always able to find a value of k/T to satisfy eqn. pV.24p for any given positive lepton asymmetry 
For 1^1 10^2 consistent with the WMAP and BBN data|63], we obtain 



T 



3 X 10"^ C - 3 X 10" 



(IV.25) 



Note that the asymmetry term from (sm) contribution dominates over the (bsm) contribution, which is different from 
/i = case where (bsm) contribution would dominate as shown in pV.17|) . This analysis leads us to conclude that for 
a lepton asymmetry hidden in the neutrino sector compatible with the bounds from ref. 62] there is the possibility of 
two MSW resonances. 



V. IMAGINARY PARTS: WIDTHS FROM VECTOR AND SCALAR BOSON DECAY. 

The quasiparticlc widths ri_2(fc) are given by eqns. pil.57[) . pil.60p . Analyzing the explicit expressions for the 
imaginary parts of the (sm) and (bsm) contributions given in the appendix, equations IA.HIA.3l and IB.HIB.3l respec- 
tively, the "on-shell" contributions are obtained from those whose 5 function constraints can be satisfied for cj ~ fc. 
It is straightforward to find that only the terms with 5{uj + p — W^p-+fc) have non-vanishing support for oj ~ k. These 
terms are given in the last lines of IA.2I and IA.3I for (sm) and the last lines of IB. 21 and IB. 31 for (bsm) contributions. 

These contributions to the quasiparticlc widths in the medium arise from the decay of the intermediate boson, either 
the vector bosons in the (sm) contributions or the scalars in the (bsm) contributions. This is depicted in Fig. pH)) . 
the Cutkosky cut through the intermediate boson (vector or scalar) yields the imaginary part. The process that 
contributes on shell a; ~ fc is the decay of the boson into the fermions (neutrinos and or charged leptons) depicted in 
this figure. 

The fact that the decay of a heavy intermediate state leads to a width was recognized in ref. [s^ . 
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k 



FIG. 10: The Cutkosky cut for imaginary part of the (sm) and (bsm) contributions, and the contribution on the mass shell 
w ~ fc. 



The analysis of the different cases is simplified by introducing 
raa(fco,A;)= Im 

Tss{ko,k)= Im 
in terms of which (see eqn. piL32|) ') 



2k 

(fcp + hk) 
2k 



Sj = 2k 



^aaiko,k) +Tssiko,k) 



A/(fco,fc) 



2k 



(V.l) 
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(V.3) 



We need these quantities evaluated on the "mass shell", namely for positive energy ko — uj{k) ^ k + M /2k. We 
find: 



h = 1 



h = 1 



Taa{k) ~ Im 



^,,(fc)~Im(AJ^-Bi^) 



(V.4) 



(V.5) 



In the above expressions we have used Yi^2 <C g and M /4fc^ <C 1 and neglected terms accordingly, we have 
suppressed the arguments on A, B, however, these matrix elements depend on k. The term with A^ + B^ in (IV. Sp is 
noteworthy: the leading contribution to this term is from (sm) interactions, even setting the Yukawa couplings in the 
(bsm) sector to zero a nearly right-handed sterile neutrino is produced via the decay of the vector bosons. 

The expression for the imaginary parts (jIII.57IIII.60[) simplify in two relevant limits [6^: 

a) weak damping: I7I <C 1: in this limit we find 

r sin((^) ~ 7 cos 29m, 
leading to the following results for the poles with positive energy 

ri(fc) - T aa{k) cos^ 9,n + sin^ e„,r ssik) , 



T2{k) = Taaik)smHm + cos^Omr,,{k). 
Furthermore the difference in the dispersion relations becomes 



(V.6) 

(V.7) 
(V.8) 

(V.9) 



22 



which is the usual result for neutrino mixing. 

b) strong damping: \^\ ^ 1, in this limit we find 



leading to the following results 



rsin(0) ~ 7 sign(cos(26'„i)) 
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(V.IO) 



ri(fc) = ^ (raa(fc) + Tss{k) ) + ^{Taaik) - Tssik) ) sign(cos(2^^„0 ) - ^'""^^^""^ 



(V.ll) 



r2(fc) = -(r,,(fc) + r,,(fc): 



liTUk) - Tssik) ) (^sign(cos(2^^„0) - 



In this case the frequency difference between the propagating states becomes 



Ailsd = n2{k) ~ n,{k) 
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(V.12) 
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This is a remarkable result, the frequency difference vanishes at a MSW resonance in striking contrast with the usual 
quantum mechanics description of neutrino mixing and oscillations wherein there is a "level repulsion" at an MSW 
resonance that prevents level crossing. 

In all the expressions above Vaa{k)] Tss{k) are given by (IV. 4[ IV. SP in the respective cases h — =Fl, and the mixing 
angle 9m is obtained from eans. (jIII.47p evaluating A/j at fco = k. 

The widths for negative energy and h = =pl are obtained from the expressions above by the replacement /x — /x, 
this is a consequence of the relations (jIIL61|IIL62|l and the fact that the chemical potential is CP-odd, therefore the 
particle and antiparticle widths only differ because of the chemical potential. 

We emphasize that the results (|V.7IV.8|) and (|V.lllV.12p are general, and hold to all orders in perturbation theory 
as they follow from the general form of the self-energies. In particular these relations are valid beyond the one-loop 
order studied here and hold for any processes that contributes to the absorptive parts of the self-energy at one-loop 
or higher order. 



A. Widths from scalar and vector boson decay: 



As discussed above, the imaginary parts of the self-energy are given in the appendix, both for (sm) and (bsm) 
contributions. Inspection of the different delta functions shows that the only contribution "on-shell" , namely uj k 
arises from the terms with S{uj + p ~ W_^_^j^) in the expressions for the imaginary parts (|B.2IB.3p . 

This delta function corresponds to a Cutkosky cut that describes the process of a scalar (in (bsm)) or a vector (in 
(sm)) boson decay into a neutrino and another lepton, displayed in Fig. (IIC 



Scalars (bsm): For scalars the (R) and (L) components are the same. We find for Y = Yi ^; Af — Ma-.tp for a,ip 
exchange respectively 



lm{Aii - Br) = Im(Ai 
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where 
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for E^^^ 



(V.16) 



23 



In the relevant region k < Ma-.ip '-^ T we can safely neglect the contribution from the leptonic chemical potential in 
(|V.14p and set ^ = 0, since the bounds from ref.jG^I suggest that |^| < 0.02. The result (|V.14p agrees with that found 
in ref.[43] for the decay of the scalar boson into sterile neutrinos (2) (Ci = 0) for vanishing chemical potential. 

For /c/T <C M/T ^ 1 we can approximate 



y2 ^2 

Im(Afl - Br) = Im(Ai - Bl) = e 



(V.17) 



Vector bosons (sm): For (sm) vector boson exchange (only L), the imaginary parts "on-shell" are obtained from 
the terms with 6{uj + p — W^p-_(_fe) in the imaginary parts (|A.2IA.3p setting ui ~ k. We find 



Im(Ai-Bi) = %^-^ln 
Idtt k^ 



1 - e-'^'-y 



(V.18) 



where gsm is given by eqn. piI.16p . M = Mz,w for neutral and charged current contributions respectively, and 



^ J 1 for SncjScc , 
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For positive energy, and positive helicity (right-handed), we also need (see eqn. (jV.sp ) 
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where Li2 is the dilogarithm or Spence's function and we have set fj, ^ 0. This expression simplifies in the limit 
fc/T < Af/T - 1 with the result 
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In the above results for vector bosons M = Mw,z respectively. 

We can now gather all the results needed for Taa{k); ^ss{k) (|V.4IV.5p and the quasiparticle widths ri.2(fc) obtained 
from them. Approximating cos(0) ^ 1; sin(0) ^ we find. 
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In the expressions above we have defined 

hp 

< = ; a = a,^,Z,Ty. (V.26) 

For smaU values of the arguments Li2{z) ^ z which may be used appropriately whenever x* > 1, a situation which 
describes the relevant range Ma T;k < T. 

Eqns (|V.22IIV.25|) combined with (|V.4IV.5|) yield the complete expressions for the quasiparticlc widths ri,2 in all 
cases, and as per the discussion below, the production rates. 



B. Imaginary parts: from the width to the production rates. 



The connection between the quasiparticlc widths (imaginary part of the self-energy "on-shell" ) and the production 
rate is established via the Boltzmann equation for the production of a given species, in this case that of a "sterile" 
neutrino. Consider as an example the scalar vertex YiVsai'a, the analysis is similar for the other, including (sm) 
vertices. The Boltzmann equation is of the form (gain)— (loss) (see for example the appendix in ref. [43|)). The gain 
term corresponds to the decay process a ^ Va + I's and is of the formjisj 
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where NbtUp are the bosonic and fermionic distribution functions respectively. The loss term describes the inverse 
process, namely the recombination Va -\- I's ^ with 
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Therefore the Boltzmann equation is of the form 
dns{k) f d^p 



dt 
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5{W.^^^-p-k)lNB{W.^j:){l-nF{p)){l-ns{k))~ 1+Nb{W.^j:) nF{p)n,{k)\ . (V.29) 



If the distribution function of the particle in question is slightly perturbed off equilibrium, the relaxation rate of 
the distribution function towards equilibrium is obtained by writing ns{k) = rfj^{k) + 5ns{k) and linearizing the 
Boltzmann equation in ^^^(fc) 59]. The linearized Boltzmann equation reads 



dSns{k) 
Jt 



-Trel Sns{k) , 



where 



I^re/ — 



(27r)3 



M 



S{W^-^i:-p-k) nFip)+NB{W^~^z) 



(V.30) 



(V.31) 



As discussed in ref. [59|, the relaxation rate Trei is twice the quasiparticle width [59| since the distribution function is 
bilinear in the fields. The relation between r^e; and the on-shell width becomes evident comparing the expression 
(IV.3ip with the "on-shell" imaginary parts, namely the last lines in eqns. (|A.2IA.3IB.2IB.3p with uj ~ k. The 
production rate of the sterile species is obtained by neglecting the inverse process and neglecting the sterile population 
buildup in the Boltzmann equation (jV.29|l . namely. 
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Therefore by obtaining the bosonic and fermionic contributions to the quasiparticle widths as in the previous section, 
we can obtain the production rate. Although the term with the product Nsnp is not included in the width, such 
term is smaller than the term with Nb only, since p'^np{p) features a maximum at p/T 2.3 for which np(p) ~ 0.09 
(for /i/T <C 1). Therefore in the region of importance in the integral p ^T, the production and relaxation rates only 
differ by a few percent, and the results for the relaxation rates yield a reliable approximation to the production rate. 
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An important bonus of obtaining the production rate from the quasiparticle decay width as carried out here is the 
correct dependence on the mixing angle in the medium, which would be missed by a naive perturbative calculation. 

Therefore the quasiparticle width yields an excellent approximation to the production rate, in particular it describes 
correctly the dependence on the mixing angles in the mediurn, its magnitude and fc-dependence. 

In particular, the result (jV.23|) confirms the result of ref.[4j] for Yi = 0. For the scalar contribution (bsm) the 
right and left-handed yield the same result (multiplying (IV.23P by a factor 2 in the total rate) and as discussed above 
the production rate is twice the width, which restores the factor 4 between (jV.23|) and the result in ref. [isj which 
corresponds to the case Yi = 0. 

Thus we conclude that the results of eqns. (|V.4IV.5IV.7IV.8IV.11IVTT2)) along with the explicit forms (|V.22IIV.25|) 
provide a complete and reliable assessment of the production rates ready to be input in the kinetic equations that 
include the cosmological expansion [43|. 



C. Weak or strong damping? 



We have now all the ingredients to assess under which circumstances the weak (I7I <C 1) or strong (I7I ^ 1) damping 
conditions are fulfilled. In terms of the widths and real parts it follows that 
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For h = 
we find 



-1, Taa — Fss and Aji are dominated by the (sm) contributions, therefore from eqns. (|IV.13|) and (|V.24p 
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where A{k),B{k) can be read off (jIV.13IV.24p . In the region of parameters where Aii{k) ^ cos(26') ^ 1, it follows 
that 7 ~ A/(fc)/Afl(A:), furthermore, for A: < T Mzw the function B{k) ^ e"^^' < 1, leading to Aj/An < 1 
corresponding to the weak damping case in which the widths (production rates) are given by (jV.7IV.8p . 

Far away from the MSW resonances but in the region where cos(20) ~ 1 ^ Aji{k) it also follows that Aj/A/j <C 1, 
corresponding again to the weak damping regime. Therefore the parameter region far away from MSW resonances 
(either above or below) corresponds to the weak damping regime. 

Very near MSW resonances cos(20) + A^ ~ and 7 ~ A//| sin(26')|, in the region of relevance for our analysis 
T Mzw with Ms ^ KeV it follows that 
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therefore, since the resonance occurs at k/T < 1 for Mz^w ~ T we conclude that the strong damping condition 7^1 
is fulfilled near MSW resonances. Because the MSW resonance(s) are very narrow for T ~ Mz^w discussed above 
(see the discussion leading to eqn. IIV.16P , we conclude that in most of the regime of temperatures and momenta the 
weak damping results (jV.TIV.Sp are valid and only in a very narrow region near MSW resonances the strong damping 
results (IV.lllV.i2p are valid. 

An identical analysis confirms a similar conclusion for the case h ^ I, namely the weak damping condition holds 
in most of the relevant range of M/T; k/T but for a narrow region near the MSW resonances in which the strong 
damping condition holds. 

An alternative interpretation of the weak and strong damping regime is obtained using eqn. (|V.9p to write 
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Since Aflsd < An^d the denominator gives an upper bound to the oscillation frequency between the active and 
sterile neutrinos. The weak damping regime I7I <IC 1 describes the case in which there are many oscillations before 
the overlap amplitude is suppressed, whereas the strong damping regime describes the case in which damping occurs 
before oscillations take place. For a similar discussion see the second ref. in[5^. 
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D. Regime of validity of perturbation theory. 

In the relativistic approximation the vahdity of the perturbative expansion requires that k 3> S&sm ; "^sm ■ Since the 
weak interaction couphng constant gsm is much larger than Yi 2 we focus on the standard model contributions. 

From the expression pV.Sp and the results displayed in fig. ([7]) we see that for M^z/T > 1, it follows that 
ReSsm oc a^T since the coefficient functions A,B < 12. Therefore perturbation theory is valid for k 3> r/30, 
hence for Mw/T > 1 the resonances in absence of lepton asymmetry at 0.2 < k/T ^ 1 for 1 < Mw/T < 3 are 
comfortably within the regime of vali dity of the perturbative expansion. The lepton-asymmetry induced resonance 
ior k/T ^ Mw/T is the usual resonance and for T <C Mw the expressions (jIV.9IIIV.ll[) reduce to the results available 
in the literature [H,!!!]. In the regime k <^T < Mw the on-shell self-energies are linear in k. We see that, for ^ 0.4, 
the terms proportional to k are ^ 1 for Mw ^ 2r, hence perturbation theory is reliable within the regime of interest 
in this article. The imaginary parts are always perturbatively small because of the exponential suppression factors 

g-M^/fcT 

Perturbation theory breaks down for M <T for the small k/T region and requires a hard thermal loop resummation 
program[6l| akin to the one presented in ref. [57| in the standard model without mixing. This is well known in gauge 
theories where the gauge bosons are nearly massless on the scale t[61[. Such program is well beyond the realm of 
this study, however for M/T > 1 — 3 our results are reliable for k/T ^ aw as analyzed above. For example for 
the case M/T ~ 1 although the peak in the coefficients A,B in the self-energy occur for k/T ^ 0.07 which is not 
too large compared to a^, ~ 0.03, the position of the resonance at k/T ^ 0.2 is well within the regime of validity 
of the perturbative expansion. The validity of perturbation theory improves dramatically for M/T > 1 even in the 
low momentum region as discussed above. Therefore, we conclude that for M/T > 1 the perturbative results are 
reliable for k/T > a^, in particular the new resonances are well within the regime of validity of the perturbative 
expansion. The results for the production rates are always perturbatively small and reliable because of the exponential 
suppression factor. 



VI. DISCUSSION 



Our goal is to study the production of sterile neutrinos in cosmology near the electroweak scale when the universe is 
radiation dominated. To include the effects of cosmological expansion in the production rates and mixing angles, one 
must first replace the momentum k kp{t) = k/a{t) and temperature T — s- T{t) — TiUi/ait) where k is the comoving 
momentum, a{t) the scale factor and T^; correspond to the initial temperature and scale factor at which the kinetic 
equations are initialized. Whereas the ratio kp{t)/T{t) — k/{Tiai) is constant M/T{t) = Ma{t)/ {TiUi) grows during 
the expansion. Consider setting initial conditions at Ti < Mw, so that M/Ti ^ 1, the analysis of section pvp shows 
that there exists at least one very narrow MSW resonance even for nearly right-handed sterile neutrinos (two if a 
lepton asymmetry in the neutrino sector is included) at a value {kp{t) /T{t)) ^ < 1. For [kp{t) /T{t)') < {kp{t) /T[t)) ^ 
the analysis shows that A/j ^ 1 and 



— «0, (VI. 1) 



therefore for {kp{t)/T{t)) < {kp{t)/T{t)) ^ wc find 



(^)r... (VI.2) 



For these values of kp{t)/T{t) the mode "1" is active-like and it is produced with a weak interaction rate, whereas 
the mode "2" is sterile like and is produced with the rate similar to that of ref. [l^ plus small corrections from 
standard model interaction rates suppressed by the mixing angle in the medium ~ 0/ ^r ■ On the other hand for 
{kp{t)/T{t)) > {kp{t)/T{t))^ we found above that <C —1 leading to 0^ ~ 7i'/2, namely the mode "1" is sterile-like 
and the mode "2" is active like, with the production rates 

2 

^ aa ^ ^ ss J 



T2 ~ r,,. (VI.3) 

As the cosmological expansion proceeds eventually M/T{t) ^ 1 and the resonances disappear (in absence of lepton 
asymmetry the MSW resonances for kp{t)/T(t) < 1 disappear for M/T{t) > 3), A^j remains large but positive and 
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the mixing angle in the medium is given by (jVI.ip and the production rates are given by (|VI.2p for all values of 
kp(t)/T{t), namely the mode "1" remains the active-like and the mode "2" the sterile-like. 
We note that (see section ()V A[) ) 



/r2 



rss,rQQ cx ( ) In 



(VI.4) 



this is precisely the form of the production rate that leads to a distribution function after freeze-out that is enhanced 
at small momentum, a feature that leads to a larger free streaming length and transfer function at small scales [43]. 

During the time when M/T{t) ~ 1 the MSW resonance for kp{t)/T{t) < 1 leads to a non-thermal population of 
neutrinos: for {kp{t) /T{t)) < [kp{t) /T{t)) ^ there is a large production of mode "1" leading to large populations and 
a small production of "2" (sterile like) leading to small populations, whereas for (^kp{t) /T{t)) > {kp{t) / T [t)^ ^ there 
is a "population inversion" in the sense that mode "1" is slightly populated whereas mode "2" will be substantially 
populated, however, without the small momentum enhancement. Consider a fixed value of kp{t) /T{t) < 1 during the 
cosmological expansion the ratio M/T{t) cx a{t) increases sweeping through the MSW resonance, when this happens 
the mixing angle in the medium vanishes very rapidly because the resonance is very narrow and the mode "2" becomes 
sterile like. As the expansion continues the MSW resonances (in absence of lepton asymmetry) disappear altogether 
and the mixing angles and production rates are given by (|VI.llVI.2jl respectively for all values of kp{t)/T{t). The 
population of the active- like neutrino (mode "1") continues to build up via weak interaction processes, including those 
that become dominant at T ^ Mw and eventually thermalizes, whereas the population of the sterile- like neutrino 
will be frozen-out as the production rate r2 shuts-off as Tss vanishes rapidly for M/T(t) <^ 1 (see ref.ji^) and 9m ^ 
as M/T{t) ^ 1 even when Taa (weak interaction rates) remain large down to the decoupling temperature of weak 
interactions ^ 1 MeV. 

This analysis indicates that sterile neutrino production via the decay of scalar or vector bosons will be effective 
only in a region for M\y /T{t) ^ 1 and the distribution function at freeze-out will be strongly non-thermal with very 
small population but with an enhancement at small momentum as found in ref. [l^. However, the weak interaction 
contribution will freeze out much later, depending on the temperature dependence of the mixing angle in the medium 
and will eventually merge with the non-resonant (DW) production mechanismjl^l at T ~ ISOMeV. 

However the non-thermal distribution built up during the stage when scalar and vector boson decay dominate the 
production will remain. 

At this stage it is important to understand the self-consistency of the analysis. In obtaining the self-energies we 
had assumed that the eigenstate "1" is active like with a thermal distribution function. We have learned, however, 
that there are resonances and the eigenstates "1" and "2" are either active-like or sterile-like depending on fc, namely 
on which side of the MSW resonance the wavevector lies. This finding calls into question the thermal nature of the 
neutrino propagator in the intermediate state (of course there is no such ambiguity in the charged lepton propagator 
that enters in the charged current self-energy). This issue notwithstanding, we have found that the fermionic and 
bosonic contributions to the real parts of the self-energies are qualitatively the same with a rather small quantitative 
difference, both for (sm) and (bsm) contributions. Therefore replacing the thermal fermion propagator for a vacuum 
one leads to a minor quantitative modification of our arguments. However because of the enormous pre-factors the 
conclusions about the sharpness of the resonance and the resonance positions do not change and the general analysis 
remains the same. Therefore, we conclude that the results obtained above are very robust not depending on whether 
the intermediate fermion line features a thermal or vacuum propagator or non-thermal propagator interpolating 
between these two cases. 



VII. CONCLUSIONS AND COSMOLOGICAL CONSEQUENCES 

A comprehensive program to assess the viability of any potential (DM) candidate begins with the microphysics of 
the production and freeze-out process of the particle physics candidate. This initial step determines the distribution 
function at freeze-out which in turn determines, along with the mass, its abundance, free streaming length, phase 
space density at decoupling and the transfer function and power spectrum in the linear regime. Our objective is to 
carry out this program for sterile neutrinos with mass in the KeV range which seems to be the range favored not only 
as a (DM) candidate but also provide potential solutions to a host of astrophysical problems [Slj. 

In this article we focus on the first step of the program and study the production of sterile neutrinos in a temperature 
regime near the electroweak scale in an extension beyond the standard model in which the see-saw mass matrix 
emerges from expectation values of Higgs-like scalars with masses of the order of the electroweak scale. This simple 
and compelling extension which features only one scale yields rich phenomenologv [sol [Sll [3^ . The main observation 
in this article is that in this temperature range sterile neutrinos are produced by the decay not only of the Higgs-like 
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scalar as explored in refs.[33,|43| but also of the charged and neutral vector bosons of the standard model. We consider 
active and sterile species to be Dirac fermions to allow the possibility to include a lepton asymmetry hidden in the 
(active) neutrino sector consistent with recent bounds from WMAP and BBN[63|. 

The assessment of the contribution from standard model vertices to sterile neutrino production requires an analysis 
of the mixing angles in the medium and production rates. We obtain both from the study of the full equation of motion 
of the active and sterile neutrinos that input the self-energies in the medium. The real part of the self-energy (index 
of refraction) determines the dispersion relations and mixing angles in the medium, and the imaginary (absorptive) 
part determines the production rates. 

We provide a detailed analysis of the contributions from "beyond the standard model" and standard model inter- 
actions to the mixing angles, dispersion relations and production rates, thereby facilitating the analysis of different 
situations. The study of the "index of refraction" in the temperature regime near the electroweak scale has not been 
performed before and yields a wealth of remarkable phenomena. 

Our study reveals the presence of narrow MSW resonances even in the absence of a lepton asymmetry, in the 
temperature regime T > Mw for k/T < 1. For vanishing lepton asymmetry the resonance occurs at a value {k/T)c 
that depends on the ratio Mw /T with 0.15 < (fc/r)c < 1 for 1 < {Mw /T) < 3. The position of the resonance {k/T)c 
increases with Mw/T, the resonance eventually disappear for Mw S> T recovering the result valid in the Fermi limit 
of the weak interactions [H, [sj. 

Including the possibility of a (small) lepton asymmetry in the neutrino sector with a value compatible with the 
bounds from WMAP and BBN[62] yields two narrow MSW resonances in these regions, with the resonance associated 
with the lepton asymmetry occurring at fc < ^ T where /i is the chemical potential for the active species that 
determines the lepton asymmetry. 

A remarkable aspect of these results is that near these resonances the contribution of the imaginary part of the 
self-energies leads to a strong damping regime, and the difference in the propagating frequencies vanishes exactly at 
the position of the resonance, with a concomitant breakdown of adiabaticity. For Mw 3> T the MSW resonances 
that are independent of the lepton asymmetry disappear leaving only the low energy resonances associated with the 
lepton asymmetry. 

Furthermore, we have found that it is quite possible that the region of parameters of the extension (bsm) allow for 
MSW resonance for positive energy, positive helicity, namely nearly right-handed states both with and without lepton 
asymmetry. We also find that the decay of the Z", vector bosons leads to the production of nearly right-handed 
sterile-like neutrinos. 

Because the resonances are very narrow, we obtain a simple expression for the production rates (see section IVip 
that is valid in a wide range of temperatures and clearly displays the contribution from standard model and beyond 
standard model interactions. 

We have argued that in the early universe the cosmological expansion leads to a highly non-thermal distribution 
function for sterile neutrinos with an enhancement of the low momentum region k < T both as a consequence of 
the MSW resonances and the vanishing of the mixing angle and production rates as the temperature falls well below 
the electroweak scale. Furthermore, we expect that because the MSW resonances are very narrow, the cosmological 
expansion will lead to sterile neutrino production resulting in a highly non-thermal distribution with low momentum 
enhancement. The form of the production rates via scalar and vector boson decay are similar to that in ref . [isj . 
which leads us to conjecture that the distribution function after freeze-out will be enhanced in the low momentum 
region, leading to a smaller free streaming length and larger power spectrum at small scales as compared to the (DW) 
mechanism (23 • 

The next step of the program will input these results into the kinetic equations that describe the production and 
freeze-out of the sterile species from which the distribution function at decoupling is obtained. We expect to report 
on these studies in a forthcoming article. 

An important remaining question is the extrapolation of these results to T ^ Mw,z- At temperatures above the 
electroweak symmetry breaking scale the SU{2) x U (1) symmetry is restored and the vector bosons become massless at 
tree level, therefore the production channel described here shuts off. However, vector bosons acquire electric screening 
masses of order gT^Gl] and scalar bosons may also acquire thermal mass corrections of 0{Yi^2T). 

Furthermore the see-saw mass matrix also vanishes at tree level if all the mass terms arise from the expectation value 
of the Higgs-like scalar field. This high temperature regime requires a deeper understanding of radiative corrections 
to the propagators of the vector bosons, in particular the hard-thermal loop corrections 61 1. 

Understanding the possibility of sterile neutrino production in this high temperature regime entails a non- 
perturbative resummation program also for neutrinos, akin to the study in ref. [57| . This program although clearly 
interesting in its own right is far beyond the realm of our goals here and deserves a deeper study. 
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APPENDIX A: VECTOR BOSON EXCHANGE (SM) 

The (SM) self-energy contributions with the exchange of a vector boson are given by the spectral representation 
(UlLSl) with the imaginary part given by eqn. piI.14p which is of the form 



2 



7° ^lm{^,P,k) - ^ ■kl{l^{uj,p,k) 



Neglecting the mass of the neutrinos and charged leptons we find 



r 1 r 
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APPENDIX B: SCALAR EXCHANGE (BSM) 

For scalar boson exchange we find 



Iml]fcs™(tj, k) = 



(Pp 



4 J {2nYW.^~, 



l^^bsnii.^,P-,k) - 1 - k (k-pj Ill^^{uj,p,k) 



(B.l) 



where 



Ul^{cj,p, k) = [l - nAp) + NB{W.^i)\ S{u;~p~ W.^f-) + [l - nAp) + Nb{W.^z)\ ^(^ + P + 

+ \np{p) + Nb{W.^j:)\ 6{u:^p + W^^j^) + \nAp) + <5(c^ + p - W.^^) , (B.2) 



nL™(^,P, fc) - [l - nAp) + Nb{W.^-,)] 5{lo~p~ W.^^) - [l - nAp) + A^s ^(^ + P + 



(p) + 7Vs(W^,~+,-) 5(^-p + iy-^^)- nAp)+NB{W.^i) 5{u,+p-W.^^) . (B.3) 
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